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We discuss two aspects of f(R) theories of gravity in metric formalism. We first study the reasons 
to introduce a scalar-tensor representation for these theories and the behavior of this representation 
in the limit to General Relativity, f(R) — > R. We find that the scalar-tensor representation is well 
behaved even in this limit. Then we work out the exact equations for spherically symmetric sources 
using the original f(R) representation, solve the linearized equations, and compare our results with 
recent calculations of the literature. We observe that the linearized solutions are strongly affected 
by the cosmic evolution, which makes very unlikely that the cosmic speedup be due to f(R) models 
with correcting terms relevant at low curvatures. 
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I. INTRODUCTION 

In the last few years, modified theories of gravity 
of the f(R) type, where R is the Ricci scalar, have 
received much attention. These theories have the 
ability to generate late-time cosmic acceleration and 
also early-time inflation depending on the details of the 
function f{R) considered. In this sense, they represent 
an interesting alternative to the combination of General 
Relativity (GR) plus dark energy models [H, 01, and 
to other explanations for the cosmic speedup [1,111]. A 
given gravity lagrangian f(R) can, in addition, lead to 
two completely different theories depending on the vari- 
ational principle used to derive the equations of motion. 
Here we will be dealing with the most common choice, 
the metric variational formalism. The other choice, 
Palatini formalism, assumes that the connection is a field 
independent of the metric and will not be considered here 
(see [1 i 0, I, 1, HE El EE EE M EE EE 113 for details) . 

The dominant underlying philosophy in the con- 
struction of f(R) models has been that non-linear 
curvature terms that grow at low curvatures could have 
a negligible effect in stellar systems or galactic scales, 
where the curvature is assumed to be relatively high, 
but a non-trivial one in the cosmic regime, where the 
curvatures involved are orders of magnitude smaller. 
This idea has led to many different proposals, am ong 
which the Carroll et al. model f(R) = R - fi 4 /R pf 
is perhaps the most known (see also [l9j for earlier 
models). A rich debate in various arenas has taken 
place with some authors focusing on the cosmological 
aspects 1 , and others on the predictions of these theories 
at smaller scales (Newtonian and post-Newtonian 
regimes). Among the latter, some claim that models 
of the 1/R type are ruled out by solar system and/or 
laboratory experiments [IE HE HH HE HE H3 (see also 
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[28| for a discussion of R n models), and others defend 
the opposite [IE HE U EE HE H| (see also [H for a 
critical discussion of these two positions). An important 
element in this discussion was the identification of f (R) 
theories with a class of scalar-tensor theories [IE HE H3 , 
which allowed one to reinterpret the equations of motion 
and made more accessible the computation of the 
Newtonian and post-Newtonian limits. More recently, 
however, this scalar-tensor approach has been criticized 
and, apparently, shown to break down for theories close 
to GR [13|. It was then claimed that solar system 
tests do not actually rule out theories of the type 
f(R) = R + eg(R), with e a small parameter, because the 
corresponding scalar-tensor theory is not well defined in 
the limit e — > 0. This conclusion, on the other hand, 
has also been criticized recently in [IE HE H3j> where 
spherically symmetric weak field configurations have 
been studied using the original f(R) representation. In 
those works, it has been found that the post-Newtonian 
parameter 7 turns out to be 7 = 1/2, in agreement with 
the results found using the scalar-tensor approach in Q • 
Solar system experiments again seem to rule out f(R) 
theories. However, there is something unclear in those 
works that should be reanalyzed before accepting their 
conclusions. In [2E HE H3 the wrong coordinate system 
was used to carry out the calculations and identify the 
post-Newtonian parameter 7. In fact, for the discussion 
of solar system experiments and the identification of 
post-Newtonian parameters, one should use isotropic 
coordinates instead of Schwarzschild coordinates (see 
[3^ . page 1097 for details). This point and others will 
be carefully addressed here. 

In this work we will discuss in some detail the 
scalar-tensor representation of f(R) models of the type 
f(R) = R + eg(R) and their limit to GR. Firstly, 
we will motivate the introduction of the scalar-tensor 
representation and will then show that it is well defined 
even in the limit e —> 0. Secondly, we will consider the 
weak field limit for spherically symmetric sources using 
the original f(R) representation. We will show that 
the same results found in j24| using the scalar-tensor 
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approach are found here using the f(R) form. For 
our calculations we use isotropic coordinates. In the 
appendix, we also solve the linearized equations using 
Schwarzschild coordinates and transform the solution to 
isotropic coordinates. This will serve as a consistency 
check for our calculations. In order to provide a complete 
description of the dynamics of f(R) theories, in our 
derivation we will take into account the interaction of 
the local system with the background cosmology. This 
will allow us to better understand the limit to General 
Relativity of f(R) theories. We will see how, driven 
by the cosmic expansion, a theory characterized by a 
low curvature scale R e goes from a General Relativistic 
phase, in which the post-Newtonian parameter 7 « 1, to 
a scalar-tensor phase, in which 7 — > 1/2. Though we will 
be mainly discussing the weak field regime, we provide 
here the basic concepts and formulas that (hopefully) 
will be used in a forthcoming work [38], in which the 
strong field regime and the contribution of higher-order 
corrections to the linearized theory [39( will be studied. 

The paper is organized as follows. We first define the 
action and derive the equations of motion for f(R) the- 
ories. Those equations are then physically interpreted, 
which motivates the introduction of the scalar-tensor rep- 
resentation. We then discuss the conditions for the equiv- 
alence between f(R) theories and scalar-tensor theories 
and study the limit f(R) — > R. The next step is to write 
the equations of motion for spherically symmetric sys- 
tems and obtain the relevant solutions for weak sources. 
We then discuss the results obtained and the interaction 
of the system with the background cosmology. Wc finish 
with a brief summary and conclusions. In the appendix 
we solve the equations using Schwarzschild coordinates 
as a consistency check for our results. 



II. THE THEORY 

The action that defines f(R) gravities has the generic 
form 

S = 2k2 / c^4;E ^/ r 5/( i^ ) + S ™[9»», VVn] (1) 

where S , m [(7,V , m] represents the matter action, which de- 
pends on the metric g^ v and the matter fields ip m , and 
k is a constant with suitable units 2 . For notational pur- 
poses, we remark that the scalar R is defined as the con- 
traction R = g^R^, where R^ is the Ricci tensor 

iV = -d A + d x T*, + r^r p pX - r^r£ A (2) 



2 We will see later how re 2 is related to the effective Newton's 
constant [see Eq, ]33H below]. 



and r^ 7 is the connection, which is defined as 

r 7 = — (<9/35a 7 + <9 7 5a/3 - dxgpj) (3) 

Variation of ([I]) leads to the following field equations for 
the metric 

f'{R)R ilv ~f{R)g llv -W^ v f'{R)+g llv Uf l {R) = k 2 T^ 

(4) 

where f'(R) = df/dR. According to (01, we see that, 
in general, the metric satisfies a system of fourth-order 
partial differential equations. The trace of (|4|) takes the 
form 

3D/'(i?) + Rf'(R) - 2f(R) = k 2 T (5) 

This equation will be useful for the physical interpreta- 
tion of the field equations. 



A. Physical interpretation 

Let us consider a generic f(R) theory not necessarily 
close to GR and rewrite ((4]) in the form 

1 k 2 1 

Rp.v^g^vR = fir R \ T v v ~ 2 wm flwffi/'vft) ~ f{ R )\ + 

+ ■j^[V tl V u f(R)-g llv nf(R)} (6) 

The right hand side of this equation can now be seen as 
the source terms for the metric. This equation, therefore, 
tells us that the metric is generated by the matter and 
by terms related to the scalar curvature. If we now 
wonder about what generates the scalar curvature, the 
answer is in (|5|). That expression says that the scalar 
curvature satisfies a second-order differential equation 
with the trace T of the energy-momentum tensor of the 
matter and other curvature terms acting as sources. We 
have thus clarified the role of the higher-order derivative 
terms present in Q. The scalar curvature is now a 
dynamical entity which helps generate the space-time 
metric and whose dynamics is determined by ([5]). 

At this point one should have noted the essential 
difference between a generic f(R) theory and GR. In GR 
the only dynamical field is the metric and its form is fully 
characterized by the matter distribution through the 
equations G^ v — k^T^. The scalar curvature is also de- 
termined by the local matter distribution but through an 
algebraic equation, namely, R = — k 2 T. In the f(R) case 
both g^u and R are dynamical fields, i.e., they are gov- 
erned by differential equations. Furthermore, the scalar 
curvature R, which can obviously be expressed in terms 
of the metric and its derivatives using now plays 
a non-trivial role in the determination of the metric itself. 
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B. Scalar- Tensor representation 

The physical interpretation given above puts forward 
the central and active role played by the scalar curva- 
ture in the field equations of f(R) theories. However, ([5]) 
suggests that the actual dynamical entity is f'(R) rather 
than R itself. This is so because, besides the metric, 
f'(R) is the only object acted on by differential opera- 
tors in the field equations. Motivated by this, we can 
introduce the following alternative notation 



<t> = f'(R) 
V{4>) = R(4>)f-f(Rm 



(7) 

(8) 



and rewrite J4} as 



1 



[V M V„0 - g^Ucj)] (9) 



Using the same notation, §5§ turns into 

3D0 + 2V{4>) — 4>-j- = k 2 T (10) 



This slight change of notation helps us identify @ and 
dTUJ) with the field equations of a Brans-Dicke theory with 
parameter lo = and non-trivial potential V(<f>), whose 
action takes the form 

S = ^2 / d 4 x^[^R - V{4>)} + S m [g^,ip m ] (11) 

In terms of this scalar-tensor representation our in- 
terpretation of the field equations of f(R) theories is 
obvious, since both the matter and the scalar field help 
generate the metric. The scalar field is a dynamical 
object influenced by the matter and by self-interactions 
according to (fTOf . 

We would like to remark that, according to the above 
derivation, the only requirement needed to express the 
field equations of an f(R) theory in the form of a Brans- 
Dicke theory is that f'(R) be invertible 3 , i.e., that R(f') 
exist. This is necessary for the construction of V(<j>). 
Therefore, the differentiability condition f"(R) ^ that 
one finds follo wing other derivations of the scalar-tensor 
representation [22l l23l . [34], [36| is just a superfluous con- 
dition generated by the particular method used 4 . 



3 The minimum requirement for a function of one variable to be 
invertible in an interval is that the function be continuous and 
one-to-one in that interval. 

4 If f"(R) vanishes at some point Ro, we may need to choose 
among different branches to extend the solutions R(f') beyond 
Rq. In other words, R(f') may not be unique in those cases. 



The inverse problem of finding the f(R) theory corre- 
sponding to a given scalar-tensor theory of the form given 
m CP also requires an invertibility condition only. In 
this case, the equations of motion lead to R = dV/d<fi. If 
this expression can be inverted to obtain <p(R) then the 
corresponding f(R) lagrangian is given by 



f(R) = R(j>{R) - V(cf>(R)) 



(12) 



We thus conclude that the condition for the equations of 
motion of an f(R) theory to be equivalent to those of 
a u> = Brans-Dicke theory is that the function f'(R) 
be invertible. Conversely, for a given to = Brans-Dicke 
theory the condition for the equivalent f(R) theory to 
exist is that the function dV \<j>) / dep = R be invertible. 



C. Limit f(R) -> R 

Now that the scalar-tensor representation of f(R) the- 
ories has been introduced and the conditions for their 
equivalence clarified, we will consider the limit f(R) go- 
ing to R to gain some insight on their dynamical prop- 
erties. We will also show that the equivalence is, by no 
means, broken in this limit. 

Let us concentrate on models of the form f(R) = R + 
eg(R) with e an adjustable small parameter. The scalar 
field is therefore identified with 0=1 + eg'(R) and be- 
comes a constant, <j) = 1, in the limit e — > 0. The relevant 
part of cj) is therefore contained in ip — g'(R). By invert- 
ing this relation we find R = R(tp) , which allows us to ex- 
press the scalar potential as V((f>) — e[R(ip)ip— g(R[<p])] = 
eV{tp). The action (JTTJ) then turns into 

1 



S=— d^x^g (1 + ev)R - e V(<p) 



2k 2 



(13) 

As expected, the limit e — > leads smoothly to the action 
of GR. Let us now look at the equation of motion for <f> 
to see what happens in that limit. In terms of <p, (fTQ|) 
can be recast as 

e [3L> + ipR{ip) - 2g(R)\ = R(tp) + k 2 T (14) 

In the limit e — > the dynamical part of this equation, 
that involving Dip, vanishes and we recover the familiar 
(algebraic) relation R = —k 2 T. This means that ip has 
completely decoupled from the theory. As a consequence 
R becomes independent of tp and we can state with 
confidence that the limit to GR is smooth. This shows 
that the equivalence between f(R) gravities and their 
related Brans-Dicke theories holds even in the limit of 
GR, since no inconsistency is found in this limit. 

To conclude this section, we mention that from (fLij) 
we can extract valuable physical information. The 
parameter e, which we now assume to be small and 
fixed, in front of the dtp term makes apparent the 
existence of two regimes. The first one corresponds to 
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a behavior very close to GR, in which R « —k 2 T is a 
good approximation. The second regime arises when the 
terms on the right hand side are small. In that situation, 
the contribution from the left hand side can no longer be 
neglected. The behaviour of R is then dominated by <p, 
i.e., the full dynamics of (fT4|) must be taken into account, 
which means that R receives contributions from T and 
from its self-interactions (or, equivalently, from tp). In 
this situation, the theory behaves as a scalar-tensor 
theory, with ip representing the scalar degree of freedom. 
We will discuss further this point later on. 



III. APPLICATIONS 

A complete description of a physical system must take 
into account not only the system but also its interac- 
tion with the environment. In this sense, any physical 
system is surrounded by the rest of the universe. The 
relation of the local system with the rest of the universe 
manifests itself in a set of boundary conditions. In our 
case, according to ([5]) and ©, the metric and the func- 
tion /' (or, equivalently, R or </>) are subject to bound- 
ary conditions, since they are dynamical fields (they are 
governed by differential equations). The boundary condi- 
tions for the metric can be trivialized by a suitable choice 
of coordinates. In other words, we can make the metric 
Minkowskian in the asymptotic re gion and fix its first 
derivatives to zero (see chapter 4 of [40( for details). The 
function f'(R), on the other hand, should tend to the 
cosmic value f'(R c ) as we move away from the local sys- 
tem. The precise value of f'(R c ) is obtained by solving 
the equations of motion for the corresponding cosmology. 
According to this, the local system will interact with the 
asymptotic (or background) cosmology via the boundary 
value f'(R c ) and its cosmic-time derivative. Since the 
cosmic time-scale is much larger than the typical time- 
scale of local systems (billions of years versus years), we 
can assume an adiabatic interaction between the local 
system and the background cosmology. We can thus ne- 
glect terms such as f'(R c ), where dot denotes derivative 
with respect to the cosmic time. 

The problem of finding solutions for the local system, 
therefore, reduces to solving ([6]) expanding about the 
Minkowski metric in the asymptotic region 5 , and ([5|) 
tending to 

3DJ'(R C ) + RJ'(R C ) - 2f(R c ) = n 2 T c (15) 

where the subscript c denotes cosmic value, far away from 
the system. In particular, if we consider a weakly grav- 
itating local system, we can take f = f' c + <p(x) and 



5 Note that the expansion about the Minkowski metric does not 
imply the existence of global Minkowskian solutions. As we will 
see, the general solutions to our problem turn out to be asymp- 
totically dc Sitter spacetimes. 



with \ip\ < \ f c \ and \ h^ v \ < 1 satisfying 
tp — > and hfu, — > in the asymptotic region. Note that 
should the local system represent a strongly gravitating 
system such as a neutron star or a black hole, the per- 
turbative expansion would not be sufficient everywhere. 
In such cases, the perturbative approach would only be 
valid in the far region. Nonetheless, the decomposition 
/' = f' c + ip(x) is still very useful because the equation 
for the local deviation tp(x) can be written as 



3Ekp + W(f c + <p)- W(f' c ) = k 2 T, 



(16) 



where T represents the trace of the local sources, we have 
defined W(f) = R(f)f - 2f(R[f}), and W(f c ) is a 
slowly changing constant within the adiabatic approxi- 
mation. In this case, ip needs not be small compared to 
f' c everywhere, only in the asymptotic regions. 



A. Spherically symmetric solutions 

The analysis of spherically symmetric solutions can be 
carried out in several ways. One of them is to write the 
line element using Schwarzschild coordinates 



ds 2 



-B{?)dt 2 



1 



■df 2 + r 2 dn 2 



(17) 



Another possibility, more suitable for the discussion of 
observable effects, is to use isotropic coordinates [37|- We 
will follow here this second option and, for completeness, 
we will discuss the other in the appendix to check the 
consistency of our calculations. We define the line ele- 
ment 

ds 2 = -A(r)e 2 ^dt 2 + J- (dr 2 + rW) , (18) 

which, assuming a perfect fluid for the sources, leads to 
the following field equations (see [4l[ for a different ap- 
proach) 
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-(19) 



^4 )■• 
2A 



(20) 



where /' = f' c + tp, and the subscripts in tp r , f' rr , Mr- 
denote derivation with respect to the radial coordinate. 
Note also that f' r = p r and f' rr = tp rr . The equation for 
<p is, according to (fTH]) and (fT8|) . 



Aip rr 



-A 



+ 4>r)<Pr 



(21) 



W(f' c + <p)-W(fi c ) , K 2 



(3P - p) 
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Equations (JT9J) , ([20]) . and (|2Tj) can be used to work out 
the metric of any spherically symmetric system regard- 
less of the intensity of the gravitational field. For weak 
sources, however, it is convenient to expand them assum- 
ing \<p\ < f' c and A = 1 - 2M(r)/r, with 2M(r)/r < 1. 
The result is 



(p rr 



-M rr (r) 

<£r_ 
fc 

2 



2 



fc 



Vc 



1 

To 



-<Pr 



fc 



P-Vc 



-<p r 



-(3P - p)+ml<p 



where we have defined 



Rcf - f(R c 



2 _ f'(Rc) — Rcf"(Rc 



3/'W 



(22) 
(23) 
(24) 

(25) 
(26) 



This expression for in 2 was first found in [24j within the 
scalar-tensor approach. It was found there that m 2 > 
is needed to have a well-behaved (non-oscillating) New- 
tonian limit. This ex pres sion and the conclusion m 2 > 
were also reached in [42j by studying the stability of de 
Sitter space. The same expression has also been found 
more recently in [2?], H^] . 

Outside of the sources, the solutions of l|22p. (|23[) and 
dlH lead to 



<p(r) = Q- e ~ m " r 
r 



A{r) = 1- 



Cn 



r 

C 2 
r 



1 - -Q- e - m A + - 
C 2 f> ) + 6 



(27) 
r 2 (28) 



Mr),:2 " = 1 ^ ( 1 + -^rf mcr ) ~ y 7 " 2 (29) 



where an integration constant has been absorbed in 
a redefinition of the time coordinate. The above so- 
lutions coincide, as expected, with those found in [24| 
for the Newtonian and post-Newtonian limits using the 
scalar-tensor representation. Comparing our solutions 
with those, we identify 



C 2 



^fc 

Ci _ 1 
fC 2 - 3 



(30) 
(31) 



where Mq = J d 3 xp(x). The line element (fT8| can be 
written as 



ds 1 



1 _2GMe_V Lr2]dt2 
r 3 



(32) 



where we have defined the effective Newton's constant 

..2 



G 



K 



1 



^fc \ 3 

and the effective post-Newtonian parameter 



7 : 



(33) 



(34) 



This completes the lowest-order solution in isotropic co- 
ordinates. The higher order corrections to this solution 
will be studied in detail elsewhere [38| (see also [3!| for 
a recent discussion of those corrections). 

Had we used the Schwarzschild line element (TIT)) , the 
result would be (see the appendix for details) 



f 



i + 5 

r 



3 

1 + m c r 



Yif 2 \ ,n- 



v 



(35) 



^r 2 ) df 2 + f 2 dn 2 



We see that the Newtonian limit, the function in front 
of dt 2 , coincides with the result obtained using isotropic 
coordinates (neglecting higher order corrections due to 
the change of coordinates). However, the first post- 
Newtonian correction is not correctly identified if we just 
take the function in front of df 2 . It is necessary to trans- 
form to isotropic coordinates using 



2 



1 



V c 
— i 
12 



(36) 



to find the correct expression. Nonetheless, in the limit 
m c r -C 1 considered in [2^, [2(| [27| and neglecting the 
V c r 2 terms, the change of coordinates is trivial (df/dr w 
1) and ([35]) agrees with ([32]) . This is, however, just an 
accident that leads to the desired result. In general, one 
should use isotropic coordinates for the correct identifi- 
cation of the post-Newtonian parameters and take into 
account the interaction with the background cosmology 
for the right interpretation of the limit to General Rela- 
tivity (see next subsection). 



B. Discussion of the linearized solutions 



We see from (13"3"|) and (|3"4")l that the parameters G 
and 7 that characterize the linearized metric depend 
on the effective mass m c (or inverse length scale m c ~ 1 )- 
Newton's constant, in addition, also depends on f' c . 
Since the value of the background cosmic curvature R c 
changes with the cosmic expansion, it follows that f' c 
and m c will also change. The variation in time of f' c 
will induce a time variation in the effective Newton's 
constant analogous to the well-known time dependence 
that arises in Brans-Dicke theories. Actually, if we use 
the scalar-tensor representation, the analogy becomes 



6 



an identity. The length scale mj 1 , characteristic of 
f(R) theories, does not appears in Brans-Dicke theories 
because in the latter the scalar potential is V(cp) = 0, in 
contrast with ([5]). 

Let us now review how the cosmic expansion proceeds 
in f(R) models with curvature terms that grow at 
low curvatures (see also [18j, [24J, |43[ for more details). 
During the matter dominated era, the cosmic expansion 
takes place as in GR. This is so because R c and T c are 
well above the characteristic scale R e of the non-linear 
terms in the lagrangian. In this situation, we find 
R c ss —k 2 T c , as discussed in the paragraph following 
equation (fT4|) . At later times, the sources become very 
diluted and R c approaches the scale R t in which the 
left hand side of (Q3J) can no longer be neglected. The 
scalar curvature then becomes a fully dynamical object 
and the theory gets into the scalar-tensor dynamical 
regime, which triggers the cosmic speedup. Therefore, 
the cosmic expansion drives the theory from a General 
Relativistic (decelerated) phase to a scalar-tensor (ac- 
celerated) phase. Let us now consider what happens in 
the solar system. During the General Relativistic phase, 
characterized by R c ^> R e , f' c — > 1, f' c ' — > 0, we find 
that m c 3> 1, which leads to e _mcr — > rapidly (very 
short scalar interaction range). During this period of 
time, we have 7 — > 1, which coincides with the General 
Relativistic result and agrees with current observations. 
At later cosmic times, when R c ~ R e , f" > 0, the 
effective mass m c becomes finite, decays with time as 
f" grows, and eventually leads to e~ m< = r « 1 over solar 
system scales (very long scalar interaction range). We 
then find 7 — > 1/2, which coincides with the value 
corresponding to the w = Brans-Dicke theory and is 
ruled out by observations. 

The change of dynamical regime in local systems can 
also be seen by looking at the value of R. To illustrate 
this point, we will use the solutions found in 24], which 
are also valid within the sources. To lowest order, we find 



R = V r . 



4tt# 



d 3 x , p(t, g) e - me |^- g '| (3?) 



\x — x'\ 



It is not difficult to show (considering a small volume cen- 
tered about \x — x'\ — ► and integrating by parts) that 
this expression becomes R = k 2 pit, x) in the limit of Gen- 
eral Relativity (R c » R e , f' c — *■ 1, f" — ► 0, m 2 — > 00, 
and V c — > 0). As the cosmic expansion takes place and 
R c approaches the scale R e of the scalar-tensor dynam- 
ical regime (/^ 1 , /" > and growing, m 2 — > finite 
and decreasing) the dependence of R on the local matter 
density is softened by the integration in (|57|) . At a given 
point x, R is then an average over nearby sources with 
weight e - m <A x ~ x \/\x — x'\ . This shows that the depen- 
dence of R on the local matter density is less and less 
important as m 2 decreases, or equivalently, as the inter- 
action range l c — m^ 1 of the scalar field grows due to the 
cosmic expansion. Note that though the magnitude of R 



is reduced at the location of the sources, with respect to 
the value n 2 p in the limit of GR, it becomes non-zero 
at farther distances from them. In fact, away from the 
sources, (|37l) can be approximated by 



R = Vc. 



m 2 K 2 M Q e 
Airf' c r 



(38) 



which has a smooth 1/r decay when m c r -C 1. Note 
that in the limit of infinite range, m c — > 0, the curvature 
R does not depend on the local sources. This shows 
explicitly that the phylosophy that motivated the design 
of f(R) models (see the introduction) was wrong, since 
R not always behaves as in GR. 

In summary, we have seen that any theory character- 
ized by non-linear terms that grow at low curvatures will 
evolve from GR during the matter dominated era to a 
scalar-tensor theory at later times. This transformation 
takes place at both large (cosmic) and short (solar 
sytem) scales. Therefore, if the cosmic speedup were due 
to the growth of non-linear corrections in the lagrangian, 
today we should be in the scalar-tensor phase, which is 
characterized by 7 — > 1/2. This fact is in conflict with 
solar system observations ("fexp = 1 + (2.1 ± 2.3) • 10~ 5 , 
pij ) and makes it very unlikely that the cosmic speedup 
be due to such corrections in the gravity lagrangian. 



IV. SUMMARY AND CONCLUSIONS 

In this work we have studied two aspects of modified 
f(R) theories of gravity, namely, their equivalence with 
w = Brans-Dicke theories and the weak field solutions 
in systems with spherical symmetry using the f(R) repre- 
sentation. We found that the equivalence between f(R) 
theories and oj — Brans-Dicke theories only requires 
invertibility of the functions f'(R) and dV/d<fi. This re- 
sult puts forward that other well-known derivations of 
the scalar-tensor representation [22|, [H, [3(| introduce an 
artificial condition, f"(R) 7^ 0, that may lead to wrong 
conclusions [34| in the limit e — > in models of the form 
f(R) = R+eg(R). We explicitly showed that such limit is 
smooth and that one recovers GR from the scalar-tensor 
theory. 

We then worked out the spherically symmetric solu- 
tions of the linearized field equations. We carried out 
the analysis idependently using two different sets of co- 
ordinates, namely isotropic and also Schwarzschild co- 
ordinates. The results are in perfect agreement when 
one transforms from one coordinate system to the other, 
which confirms the validity of our results. We have also 
made emphasis on the adiabatic interaction between lo- 
cal systems and the background cosmology via boundary 
conditions. We have seen that the parameters that char- 
acterize the linearized solutions depend on the form of the 
lagrangian and on the value of the cosmic curvature R c at 
a given cosmic time. It is then straightforward to see that 
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the GR solution is recovered in the limit R c ^> R e , where 
R e represents the characteristic low curvature scale of the 
lagrangian. When R c , driven by the cosmic expansion, 
approaches the scale R e , there is a change of dynamical 
regime and the theory becomes closer to a scalar-tensor 
theory than to GR. This change makes the predictions of 
the theory incompatible with solar system observations, 
since 7 goes from 7 ~ 1 in the GR phase to 7 — > 1/2 
in the late stages of the scalar-tensor phase. Therefore, 
since the same mechanism that predicts cosmic speedup 
is also responsible for the running of the parameter 7, it 
seems very unlikely that this type of f(R) theories may 
be responsible for the observed cosmic acceleration. 



For weak sources it is convenient to expand (|A2[> . (|A3|) . 
and (|A4]) assuming \ip\ < f' c , D(f) = 1 - 2M(f)/f, with 
2M(f)/f < 1, and B(r) = 1 - $(f), with $(f) < 1. The 
result is 



2 



— p + V c + — I fff + -Iff 
J c Jc \ ' 

2M(r) , v - , 2^ f 



(A5) 
(A6) 
(A7) 



Appendix A: SCHWARZSCHILD COORDINATES 



Using a line element of the Schwarzschild form 

1 



ds* = -B(f)dr 



—df 2 + ? 2 dtl 



D(f) 



(Al) 



the equations that follow from ([6]), assuming a perfect 
fluid, are 



D 
~2 



11 

r 



r 



D f 
~~2 



~B 



1 -D 



1 -D 



Rf - f(R) 



2/' 



- (r- + -n 

j, \ Jrr 1 r Jr 



f r 



Rf - /Cg) 
2/' 



(A2) 



(A3) 



where V c and m 2 have the same definitions as in the 
isotropic case. The solutions to the above equations lead 
to 



(p(r) 
M(r) 



(A8) 



2C 2 

f 



1 



|f 2 (A10) 



Comparing the Newtonian potential $(f) with the 
isotropic solution, we identify C2 = 2C 2 , C\ — C\. A 
transformation from Schwarzschild to isotropic coordi- 
nates [see (|3"6| ] shows that this identification is consistent 
with the first post-Newtonian correction, and confirms 
the validity of our calculations. 



where /' = f' c + tp, and the subscripts in Df, Bf, f~, f-~ 
denote derivation with respect to the radial coordinate. 
Note also that f~ — ipf and f~~ = ip ff . The equation for 
ip is, according to (Tl7j)) and (|A1|) . 
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